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The chemical abundance patterns observed in metal-poor Galactic halo stars contain the 
signature of the first supernovae, and thus allows us to probe the first stars that formed in the 
universe. We construct a theoretical model for the early chemical enrichment history of the 
Milky Way, aiming in particular at the contribution from pair- instability supernovae (PISNe). 
These are a natural consequence of current theoretical models for primordial star formation at 
the highest masses. However, no metal-poor star displaying the distinct PISN signature has yet 
been observed. We here argue that this apparent absence of any PISN signature is due to an 
observational selection effect. Whereas most surveys traditionally focus on the most metal-poor 
stars, we predict that early PISN enrichment tends to 'overshoot', reaching enrichment levels 
of [Ca/H] ~ —2.5 that would be missed by current searches. We utilize existing observational 
data to place constraints on the primordial initial mass function (IMF). The number fraction 
of PISNe in the primordial stellar population is estimated to be < 0.07, or < 40% by mass, 
assuming that metal-free stars have masses in excess of 10 Mq. We further predict, based on 
theoretical estimates for the relative number of PISNe, that the expected fraction of second 
generation stars below [Ca/H] = —2 with a dominant (i.e., > 90%) contribution from PISNe 
is merely ~ 10~ 4 — 5 x 10~ 4 . The corresponding fraction of stars formed from gas exclusively 
enriched by PISNe is a factor of ~ 4 smaller. With the advent of next generation telescopes and 
new, deeper surveys, we should be able to test these predictions. 

Subject headings: cosmology: theory — galaxies: high-redshift — stars: abundances — stars: Population 
II — Galaxy: evolution — Galaxy: halo 
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Little is known about the first stars in the uni- 
verse, the so-called Population III (Pop III hence- 
forth). Over which range of masses were they 
formed? And what was their relative distribution 
by mass, i.e., their initial mass function (IMF)? 
The apparent lack of metal-free stars in the Galac- 
tic halo suggests that star formation in metal- 
free gas should be skewed towards relatively high 
masses, as compared to the mass scale of < 1 Mq 
found in the solar neighborhood (e.g., Bond 1981; 
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Oey 2003; Karlsson 2005; Tumlinson 2006). This 
finding is corroborated both by theoretical argu- 
ments (e.g., Larson & Starrfield 1971) and by nu- 
merical simulations (Bromm et al. 1999; Bromm 
et al. 2002; Abel et al. 2002; Yoshida et al. 
2006) following the collapse of primordial bary- 
onic matter inside the very first non-linear struc- 
tures, virialized dark matter halos of ~ 10 6 Mq, 
so-called minihalos. As a consequence of the inef- 
ficiency of primordial gas to cool to temperatures 
below ~ 10 2 K, these simulations indicate that 
the first stars were predominantly very massive 
objects with masses > 100 Mq. 



It can be argued that star formation in mini- 
halos, apart from being regulated by H2 cool- 
ing, occurred under simplified conditions where 



the absence of magnetic fields, ultraviolet radia- 
tion, and a pre-established turbulent velocity field 
resulted in a unique metal-free stellar population 
(e.g., Bromm & Larson 2004; Glover 2005). For in- 
stance, in shock-compressed or photo-ionized pri- 
mordial gas, the formation of the coolant HD is 
boosted, leading to a decrease in the typical mass 
of metal-free stars down to ~ 10 A4® (Johnson 
& Bromm 2006). Similarly, lower mass stars may 
also be formed in turbulent molecular clouds with 
a negligible magnetic field (Padoan et al. 2007). 
The presence of a magnetic field could also, in 
the case of magnetized accretion disks, lead to 
smaller stellar masses as the accretion onto the 
proto-stellar core may be reduced by magnetically 
driven outflows (Silk & Langer 2006). 

In spite of their peculiarity and apparent rar- 
ity, the first, very massive stars ultimately set the 
stage for subsequent star and galaxy formation. 
If mass-loss through rotationally induced stellar 
winds was limited (Ekstrom et al. 2006), a signif- 
icant fraction of these stars was likely to explode 
as pair-instability supernovae (PISNe), leaving no 
remnant behind (Heger & Woosley 2002). Inter- 
estingly, i] Carinae and the Pistol Star, both in the 
mass range ~ 100 — 200 A4q, are two examples of 
such very massive stars in our own Galaxy. More- 
over, the recent detection of a particularly ener- 
getic supernova (SN) explosion, SN 2006gy (pro- 
genitor mass estimated to m ~ 150A^q) in NGC 
1260, indicates that stars may explode as PISNe, 
even at low redshift (Smith et al. 2007). 

Since very massive stars are extremely short- 
lived (< 3 Myr), primordial PISNe (140 < 
tu/Mq < 260, Heger & Woosley 2002) would sup- 
posedly be the first stellar objects to enrich the 
interstellar medium (ISM) with metals. As pre- 
dicted by homogeneous chemical evolution models 
(e.g., Ballero et al. 2006), such an initial en- 
richment would establish a metallicity floor and 
imprint the chemical signature of PISNe in the 
early ISM, from which the first low-mass (m < 
1 M ) stars later formed. The unique chemical 
signature of the primordial PISNe should thus 
predominantly be retained in the oldest and, in 
particular, the most metal-poor stars belonging to 
the Galactic halo population. To date, no such 
signature has been found (e.g., Christlieb et al. 
2002; Cayrel et al. 2004; Cohen et al. 2004; 
Barklem et al. 2005; Frebel et al. 2005). This has 



been regarded as an indication that stars in the 
early universe with masses in excess of 100 A4q 
were extremely rare (e.g., Tumlinson et al. 2004; 
Ballero et al. 2006), if not altogether absent. How- 
ever, we shall argue that true second generation 
stars which were formed from material predomi- 
nantly enriched by primordial PISNe have metal- 
licities significantly above those of the stars in the 
metal-poor tail of the Galactic halo metallicity 
distribution function (MDF), at variance with the 
classical picture described above. Furthermore, 
although the mass fraction of Pop III stars ending 
up as PISNe may be significant, the fraction of 
low-mass stars bearing their chemical signature is 
small, which explains the absence of such stars in 
previous surveys. 

The organization of the paper is as follows. In 
£j2] we consider the chemical enrichment by the 
first stars in its cosmological context, while the 
detailed chemical enrichment model and parame- 
ters are discussed in Sj3] The results are presented 
in 21 and we conclude with a discussion of the 
implications in Sj5j 

2. The Cosmological Context 

The end of the epoch known as the cosmic dark 
ages, when our universe witnessed the first sources 
of stellar light, is one of the final frontiers in mod- 
ern cosmology. One indirect way to search for the 
first, Pop III stars is to detect their contribution to 
the cosmic infrared background. Although tanta- 
lizing progress has been made (e.g., Kashlinsky et 
al. 2007) , we have to await the advent of next gen- 
eration telescopes, such as the James Webb Space 
Telescope (JWST) and the Atacama Large Mil- 
limeter /submillimeter Array (ALMA) (e.g., Wal- 
ter & Carilli 2007) , to be able to resolve individual 
sources and their birth sites, and to place stronger 
constraints on the IMF and the star formation rate 
(SFR) in primordial galaxies. 

A complementary approach is to search for 
Pop III signatures in the local fossil record, i.e., to 
look for specific chemical abundance ratios in the 
atmospheres of old, low-mass stars in the Milky 
Way which may have formed from gas enriched 
with the ejecta from the first SNe. This gen- 
eral strategy is termed near-field cosmology, and 
this particular approach is a powerful way to learn 
about the first generations of stars. Unlike stars 
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exploding as core collapse SNe, stars in the PISN 
mass range (i.e., 140 < ttl/AAq < 260) exhibit 
only a very small neutron excess in their interiors 
(Heger & Woosley 2002). As a consequence, there 
is a pronounced odd-even effect, i.e., particularly 
low abundance ratios of odd-Z elements to even- 
Z elements, in the ejecta of PISNe. Furthermore, 
due to the lack of excess neutrons in addition to 
less-rapid expansion timescales during the explo- 
sion, the rapid n-capture process presumably does 
not take place in PISNe (Heger & Woosley 2002; 
Umeda & Nomoto 2002). These and other chem- 
ical characteristics observable in low-metallicity 
stars, such as low values for Fe/Ca, could, in prin- 
ciple, be used to identify a possible pre-enrichment 
by PISNe. However, despite considerable efforts 
made to identify and analyze the metal-poor tail 
of the Galactic halo MDF (sec the in-depth review 
by Beers & Christlieb 2005), no observational ev- 
idence for a population of very massive stars ex- 
ploding as PISNe has yet been found. It may ap- 
pear, therefore, that only a negligible fraction of 
metal-free stars had masses in excess of 100 Mq. 
We will here argue, however, that this apparent 
absence of PISN progenitors might arise from a 
subtle observational selection effect. 

Classical, homogeneous chemical evolution 
models (see e.g., Prantzos 2005; Ballero et al. 
2006) predict that signatures of primordial PISNe, 
if they ever existed, by necessity must dominate 
the chemical abundance patterns found in the 
most metal-poor stars for two reasons: the very 
short lifetimes of the PISN progenitors and the 
assumption that SN ejecta are mixed instanta- 
neously within the entire system. Evidently, such 
calculations would suggest that all true second 
generation stars, i.e., stars enriched by a single 
primordial SN, will have the lowest metallicities 
of the observed metal-poor stars. This, however, 
may not be the case. As we will show, the chemical 
signature of primordial PISNe may only be found 
in stars of somewhat higher metallicity, possibly 
explaining why searches for the most metal-poor 
stars have not turned up any stars showing the 
PISN signature. 

Although atomic diffusion in general is boosted 
by turbulence, giving rise to turbulent diffusion, 
the mixing of mass in the ISM is not instanta- 
neous. This fact becomes increasingly important 
in lower metallicity systems, as it implies that 



there should exist coeval regions in the ISM en- 
riched by different types of SNe with distinctly 
different abundance patterns. In particular, the 
amount of, e.g., calcium produced in a PISN can 
be several orders of magnitude higher than what 
is synthesized in a normal core collapse SN, with 
ratios such as Fe/Ca that are only ~ 1/10 of those 
in core collapse SNe. A second generation star 
formed out of gas enriched by a single primordial 
PISN may therefore have a metal content which 
is atypical and significantly higher than a star 
formed from gas enriched by one, or even several 
core collapse SNe (see Bromm et al. 2003; Greif 
et al. 2007). This effect is illustrated in Figure Q] 
which shows, in a statistical sense, the spatial vari- 
ation in the level of chemical enrichment that may 
be expected in the early universe. Also shown is 
the corresponding variation in the number of SNe 
contributing to the enrichment. Notice the poor 
correlation in certain regions between this number 
and the actual level of enrichment. As a result 
of the non-instantaneous mixing of metals, which, 
by construction, is not handled in homogeneous 
chemical evolution models, the chemical relics of 
the very first stars in the universe (i.e., the pri- 
mordial PISNe) could predominantly be "hidden" 
in a population of relatively metal- rich stars. 

Since we are unable to follow the hierarchical 
build-up of the Milky Way halo in detail, we shall 
focus on the chemical enrichment of individual so- 
called 'atomic cooling halos' of mass ~ 10 8 A4q, 
which are able to cool due to lines of atomic hy- 
drogen. These halos are probably massive enough 
to harbor continuous star formation while small 
enough to be dispersed during the subsequent as- 
sembly of the Galactic halo, either one-by-one or 
by first being incorporated into larger halos (see 
Robertson et al. 2005). Likely, a significant frac- 
tion of the metal-poor stars in the Galactic halo 
originate from such atomic cooling halos, in which 
the first SN-enriched gas may have re-collapsed 
and formed Pop II stars (e.g., Greif et al. 2007), 
including the second generation stars that we seek 
to identify. In accordance with this picture, Scan- 
napieco et al. (2006) demonstrated that second 
generation stars should be expected over a wide 
range of Galactocentric radii. We shall assume 
that the star formation history of the atomic cool- 
ing halos begins with the first stars formed in 
collapsed minihalos. This picture is in agree- 
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Fig. 1. — Illustration of the inhomogeneous chemical enrichment model. SNe of different types enrich space 
during the assembly of a 10 8 7W© halo. Left panel: The number of contributing SNe in the interstellar 
medium at some early time. Right panel: The corresponding Ca abundance. Although the ejecta of PISNc 
are spread over very large volumes, the Ca abundance in regions enriched by a single PISN is significantly 
higher than in regions enriched even by numerous core collapse SNe. 



ment with hierarchical structure formation, where 
smaller systems merge to form larger ones. In the 
following, star formation will be divided into two 
modes: Pop III star formation, which only occurs 
in metal-free gas, and Pop II star formation, oc- 
curring in gas already enriched with metals. 

3. Modeling Early Chemical Enrichment 

In order to discern the chemical signature of the 
first massive Pop III stars, we have constructed 
a model for how the early universe was enriched 
with heavy elements ejected by the first SNe. Our 
model accounts, self-consistently, for the forma- 
tion of both Pop III and Pop II stars, and allows 
us to study in some detail how SNe from each of 
these populations contributed to the abundances 
that can be observed today in Galactic metal-poor 
stars. In this section, we describe our chemical 
evolution model, emphasizing how it incorporates 
the distinct characteristics of Pop III star forma- 
tion in a cosmological setting. 

3.1. The Basic Model 

Although the details of the inhomogeneous 
chemical enrichment model are described else- 



where (see Karlsson 2005, 2006; Karlsson & 
Gustafsson 2005) , we will provide a general outline 
and discuss relevant improvements below. The 
model is based on the mixing volume picture, 
where the metal-rich ejecta of individual SNe are 
spread within finite but continuously growing vol- 
umes V m i x (t). In this picture, the mixing volume 
growth is modeled as a random- walk process, pre- 
dominantly driven by turbulent gas motions. At 
any given time t, we may thus define the size of 

V^nix aS 

47T 

VmUt) = —{W + CJ E fl 2 , (1) 

where D t — (u t )/ t /3 is the turbulent diffusion co- 
efficient and ue is the minimal mixing area, intro- 
duced to account for the initial expansion of the 
ejected material, before it merges with the am- 
bient medium. Furthermore, in accordance with 
the diffusion approximation, the concentration of 
metals within V, n i X is assumed to be non-uniform, 
following a Gaussian density profile. This was not 
taken into account in earlier versions of the model. 
The turbulent diffusion coefficient is estimated to 
be D t ~ 3.3 x 10~ 4 kpc 2 Myr -1 , assuming a typ- 
ical turbulent velocity of (v t ) ~ v vil ~ 10 km 



4 



s -1 , where u v ir is the virial velocity of ~ 10 8 Mq 
atomic cooling halos, and a turbulent correlation 
length of h = 100 pc, a typical value for the virial 
radius of a minihalo, on which scale mergers may 
be expected to drive turbulent mixing. 

During the early expansion of the SN rem- 
nant itself, no star formation is expected to oc- 
cur within Vmix. This is taken into account by 
assuming an initial, instantaneous growth of V m i x 
set by o e- The value of cte depends on how much 
mass of gas is swept up by the remnant, which is 
chiefly determined by the explosion energy of the 
SN and only very weakly dependent on the den- 
sity and metallicity of the ambient medium (Cioffi 
et al. 1988). For core collapse SNe, this initial 
dilution mass is estimated to be M^w — 10 5 A4q 
assuming an explosion energy of E = 10 51 ergs, 
an average ISM metallicity of Z/Z Q — 10 -3 , and 
a density of n ~ 0.1 cm~ 3 (the weak depen- 
dences on Z and n are neglected, see Ryan ct 
al. 1996 and references therein), while for PISNe, 
M dl] ~ 5 x 10 5 (£/10 51 ergs) - 6 M & at z ~ 20 
(Bromm & Loeb 2003). Depending on the PISN 
explosion energy, the associated initial dilution 
mass may be nearly two orders of magnitude larger 
than the corresponding dilution mass for core col- 
lapse SNe. Clearly the inclusion of this effect is 
important as it may strongly influence the abun- 
dance of chemical elements in regions where sub- 
sequent stars eventually are allowed to form (see 
Fig. P). 

The present model is modified to handle the 
statistics of the metal-free Pop III and the metal- 
enriched Pop II separately. Pop II star forma- 
tion is assumed to be spatially uncorrelated. This 
means that the probability, wn(fc,t), of finding a 
region in space enriched by the ejecta of k Pop II 
SNe at time t can be approximated by the Poisson 
distribution, i.e., 

u*(fc, t) = Po(fc, ,in(i)) = e-*"<*Vn(t)Vfc!, ( 2 ) 

where /zn(t) denotes the space- averaged number 
of Pop II SNe contributing to the enrichment at 
time t. For Pop III star formation, however, the 
situation is somewhat different. If Pop III stars 
would have been randomly distributed in those 
areas where star formation occurred, i.e., in the 
higher density filaments, the probability of find- 
ing a region in such a filament enriched by a given 



number k of Pop III SNe would have been de- 
scribed by equation (J5J). However, the fact that no 
metal-free star is allowed to form in a region previ- 
ously enriched by metals (but see Wyithe & Cen 
2007; Cen & Riquelmc 2007) introduces a type 
of anti-clustering effect. The fraction of enriched 
gas will thus be larger than that described by the 
Poisson statistics, for any given value of /xni(i) > 0, 
the space-averaged number of contributing Pop III 
stars at time t (see M3.4[) . To approximately ac- 
count for this "no overlapping" -effect, the Poisson 
distribution is empirically modified by a Gaussian. 
Thus, the probability wm of finding a region en- 
riched by k Pop III stars is given by 

ww(k, t) = Cfc e-(^(*)- fe ) 2 /4 x Po(fc, m (t)), (3) 

where Ck is a normalization factor. In contrast 
to wn, the probability density function described 
by equation §3§ is more narrowly peaked towards 
/im = k. In particular, wn predicts a faster-than- 
exponential decrease of the fraction of primordial 
gas available to Pop III star formation, as a func- 
tion of /im- With time, ium(fc = 0,t) will still de- 
crease slower than predicted by the corresponding 
Poisson equation since the Pop III SN rate, in fact, 
is proportional to the available fraction of primor- 
dial g will be further discussed in < j3.4l 

Given the expressions in equation (JSJ) and ([3]), 
we may calculate the probability /m(fc' ', fc) that a 
low-mass star (here defined to be a star which 
has survived to the present, i.e., with a mass 
< 0.8 Mo) is formed out of gas enriched by a 
total number k of SNe (i.e., Pop III + Pop II), k' 
of which are Pop III SNe. This probability is given 
by the integral 

TH 

fm{k',k) = c k J a*(t)wm(k' ,t)w^(k-k' ,t)un(t)dt. 
o 

(4) 

Here, wn denotes the formation rate of Pop II stars 
expressed in terms of the resulting SN rate (see 
N3.4p . while a*(t) is the fraction of still surviving 
stars in a stellar generation formed at time t and c k 
is a normalization factor (note that c k ^ C&). The 
age of the system is set to th — 10 10 yr to account 
for late-time star formation in the Milky Way halo. 
Equation (|U) is used to predict the Pop III contri- 
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bution to the amount of metals in the gas from 
which the low-mass stars were formed. Follow- 
ing Karlsson (2005, 2006), we make use of Equa- 
tions H]) — (HJ) to calculate the predicted relative 
number density of still surviving (to < 0.8 Mq) 
metal-poor stars with any given chemical compo- 
sition. These density functions are then analyzed 
and compared to observations. In the remainder 
of this section, we will define and discuss the other 
model parameters. 

3.2. Density Evolution 

Due to the inability of our chemical enrichment 
model to account for spatial variations, except in 
a statistical sense, we will adopt the average den- 
sity evolution n — n(t) of a spherically collapsing 
10 8 Mq halo embedded in an expanding universe. 
We assume that the density of the halo follows 
the top-hat model for its collapse and virialization 
(e.g., Padmanabhan 2002; Bromm et al. 2002). 
After virialization at z v ; r ~ 10, a constant den- 
sity of n = n v i r = 0.1 cm -3 is assumed, where n v ; r 
denotes the average density at virialization. 

3.3. Initial Mass Functions and Stellar 
Yields 

The relative distribution by mass for metal- 
free and extremely metal-poor stars is largely un- 
known. This is especially true for the high-mass 
end of the IMF since all metal-free, massive stars 
in our part of the universe are long gone. We must 
therefore resort to theoretical predictions. We 
shall assume that the IMF of the metal-enriched 
Pop II follows that of the solar neighborhood, 
which can be approximated with a broken power- 
law, i.e., dN/dm = 4>{m) oc m Q , with a Salpeter- 
like exponent of a = —2.3 above to = 0.5 A4q 
and a shallower slope of a = —1.3 below this 
mass (Kroupa 2001). Very massive stars explod- 
ing as PISNe are assumed to form exclusively in 
metal-free gas and are therefore not formed in the 
Pop II mode (but see Langer et al. 2007). Con- 
versely, low-mass stars are only allowed to form in 
the Pop II mode. In addition, following Bromm 
& Loeb (2003), low-mass star formation is re- 
stricted to occur in regions with efficient cooling, 
i.e., where the abundance of carbon and/or oxygen 
is high enough. To include this effect, we adopt the 
metallicity criterion £>trans, introduced by Frebel 



et al. (2007), which takes into account the simul- 
taneous cooling by O I and C II. Hence, low-mass 
stars are only able to form in gas with C and O 
abundances such that 

Arans = log 10 (lO^ + 0.3 X lO^) > -3.5 . 

The special initial conditions in the first collaps- 
ing minihalos presumably favored the formation of 
a population of very massive stars. Although these 
conditions are not always met in metal-free gas in 
general (see, e.g., Johnson & Bromm 2006; Padoan 
et al. 2007), the top-heaviness of the primordial 
IMF seems to persist. In particular, the appar- 
ent lack of metal-free stars in the present Galactic 
halo indicates that primordial low-mass star for- 
mation was highly suppressed (e.g., Karlsson 2005; 
Tumlinson 2006). For the Pop III IMF, we shall 
assume that stars below 10 Ai® were not able to 
form. Above this limit, stars in the mass range of 
core collapse SNe are assumed to be distributed 
in the same way as for Pop II, i.e., as a power- 
law with a slope of a = —2.3. Note that although 
Padoan et al. (2007) derive a primordial/metal- 
poor IMF with a slope which asymptotically ap- 
proaches a — —3.5 at the high-mass end, the slope 
in the core collapse range is significantly shallower. 

As regards the shape of the unknown, metal- 
free very-high-mass IMF, we take a conservative 
approach and assume that stars in the mass range 
of PISNe are all formed with equal probability, 
independent of their individual masses. As a mea- 
sure of the fraction of very massive Pop III stars, 
we introduce the parameter (3. This parameter 
is defined as the ratio between the number of 
PISNe and the total number of exploding stars 
in a Pop III generation. With this definition, (3 is 
a lower limit to the original fraction of very mas- 
sive Pop III stars, as a fraction of these stars will 
become black holes by direct collapse instead of 
exploding as PISNe (e.g., Heger et al. 2003). 

In the following, j3 will be considered as a free 
model parameter, and, as we shall see, it is also 
one of the parameters to which the result is most 
sensitive. If (3 = 1, all Pop III stars were formed 
as very massive stars ending their lives as PISNe, 
while if /3 = 0, no PISN ever occurred in the early 
universe. In later sections, we shall constrain j3 by 
comparing observations with model results. It is, 
however, possible to estimate the plausible range 
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for /3 from a theoretical point of view. By extrap- 
olating the Galactic IMF (Kroupa 2001) to very 
high masses, f3 is estimated to be 0.032. Alter- 
natively, using the more realistic primordial IMF 
by Padoan et al. (2007), a /3 of 0.017 is obtained. 
Greif & Bromm (2006) estimated the mass frac- 
tion of metal-free stars that went into very massive 
stars by considering the relative number of miniha- 
los and atomic-cooling halos in the early universe. 
They derived a value of 5 — 10% by mass, which 
translates into a j3 ~ 0.005 — 0.01, based on the 
idea that the PISNe and core collapse SNe were 
formed in different types of halos, the former in 
minihalos and the latter in atomic halos. 

We consider here only those chemical elements 
which are predominantly synthesized in massive 
stars, e.g., calcium and iron. We will not take into 
account the limited enrichment by intermediate- 
mass stars and thermonuclear (Type la) SNe, 
whose main contributions to the chemical enrich- 
ment occur at a considerably later stage in the 
history of the Milky Way. The ejection of newly 
synthesized material will be restricted to two mass 
ranges: that of core collapse SNe and that of 
PISNe (Heger & Woosley 2002). The stellar yields 
of core collapse SNe, with progenitor masses in the 
range 13 < m/M.® < 40, are taken from the new 
calculations by Nomoto et al. (2006), while the 
yields of the PISNe (140 < m/M Q < 260) are 
taken from Heger & Woosley (2002). Core collapse 
SNe with progenitor masses < 13 Mq are believed 
to synthesize very few heavy elements (see, e.g., 
Nomoto 1987; Mayle & Wilson 1988) and are ne- 
glected here. As an estimate of the uncertainty 
in the results due to uncertainties in the nucle- 
osynthesis calculations (see H4.4p , we also ran a 
simulation using the yields of PISNe by Umeda & 
Nomoto (2002). 

3.4. Star Formation Rates 

As discussed above, the metal-free, Pop III 
mode of star formation likely differed from that 
of Pop II, occurring in metal enriched gas. We 
will account for this distinction by introducing one 
SFR for the Pop III mode and a different one for 
the Pop II mode. Below, these SFRs will be dis- 
played in terms of their respective SN rate in units 
of kpc -3 Myr -1 . Stars not ejecting any metals, 
e.g., stars subject to direct black hole formation, 
are not included in these rates. Again, low-mass 




Fig. 2. — Sites of Pop III star formation in our 
SPH simulation. Shown here in projection is the 
density field around the ~ 10 8 M system in our 
cosmological box at redshift z ~ 10. The green 
dots indicate the sites where the primordial gas 
collapses into minihalos and so allows the forma- 
tion of Pop III stars. We estimate the resulting 
star formation rate of Pop III stars by assuming 
that a single Pop III star forms in each of these 
minihalos. 

star formation is considered to occur in Pop II 
only, with the low-mass star formation rate closely 
following that of Pop II SNe. 

We carried out a smoothed particle hydrody- 
namics (SPH) simulation of the assembly of a 
10 8 A4q halo, including star formation without 
feedback, in order to estimate the star formation 
rate for Pop III stars forming in minihalos which 
are later incorporated into the larger halo. The 
simulation is set up similarly to those conducted 
in Johnson et al. (2007), employing a cosmologi- 
cal box of comoving length 660 kpc. We track the 
sites where the primordial gas is able to collapse 
to form Pop III stars in this simulation (Fig. [2]), 
and we thereby estimate the Pop III star forma- 
tion rate within the volume being assembled into 
the 10 8 Mq halo, assuming that a single Pop III 
star forms in each of the collapsing minihalos (e.g., 
Yoshida et al. 2006). In the simulation, the first 
star forms at a redshift of z ~ 23. We find that the 
Pop III star formation rate, here expressed as the 
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Fig. 3. — Adopted Pop III (solid lines) and Pop II 
(dashed lines) SN rates for (3 — 1, 0.1, and 
0.01 (darker lines correspond to smaller 0). The 
dashed, vertical (brown) line marks the virializa- 
tion redshift of the 10 8 yVf© halo. In general, as 
an increasing volume of space becomes metal en- 
riched, the Pop III rates decrease while the Pop II 
rates increase. The bump in the Pop III rates at 
z — 10 is due to a sharp rise in density prior to 
virialization. The small delay in the rise of the 
Pop II rates is due to the finite lifetimes of Pop II 
core collapse SNe. 

resulting SN rate, is ~ 0.04 kpc~ 3 Myr -1 (physi- 
cal) within the region which is finally virialized in 
this halo. We thus use this value for the normal- 
ization of the Pop III star formation rate. We note 
that the neglect of local radiative and mechanical 
feedback in this simulation is not likely to greatly 
affect the value that we find for the star formation 
rate (see Susa & Umemura 2006; Ahn & Shapiro 
2007; Johnson et al. 2007; Greif et al. 2007; but 
see Whalen et al. 2007). 

The star formation rate is assumed to be pro- 
portional to the gas density n(t). Hence, 



unit) = u n ,o x Q p (t) 



n(t) 



0.1 cm" 



(5) 



and 



the volume filling factor of primordial gas. The 
Pop III rate is normalized to un,o = 0.04 kpc~ 3 
Myr -1 (see above). For metal enriched gas, the 
star formation efficiency is assumed to be 10 times 
higher. The Pop II rate is therefore normalized to 
un,o = 0-4 kpc~ 3 Myr -1 which is close to the av- 
erage SN rate in the galactic halo, derived from 
chemodynamical modeling (Samland et al. 1997). 
The fraction of stars predominantly enriched by 
PISNe will not be very sensitive to these normal- 
izations (see Table [2]). We assume that star for- 
mation begins at z — 30, to account for the rare 
high-er peaks in the dark matter distribution. 

In the above expressions, the presence of the 
factor Q p is important as it regulates the amount 
of gas available for both the Pop III and the Pop II 
mode of star formation. It is controlled by the 
Pop III SFR via Equation ©. Hence, for k = 0, 



Q p (t)=wm(0,t) 



— e -AWi(t)-jum(*) 



/4 



(7) 



Since the Pop III SN rate depends on the avail- 
able amount of primordial gas which, in turn, is 
determined by the rate of metal enrichment, uu 
and Qp are coupled via the expression for /in: 



m(t) = PJ v^(t-t')u m (t')at 





+ 0—0) I V'Ut ~ t')u w {t')dt' . (8) 
o 

In this expression, the contribution from PISNe 
and core collapse SNe are distinct since cte, which 
controls the initial expansion of the ejecta, dif- 
fers between the two types of SNe (see above). 
Here, V^ m and V££ x denote the mixing volumes for 
PISNe and core collapse SNe, respectively, while 
/?, defined in §3.3| denotes the fraction of PISNe 
in Pop III. 

4. Results 



m(t) = m,o x (l - Qp(t)) 



n{t) 



0.1 cm" 3 



(6) 



where ui and denote the Pop III and Pop II 
SN rate, respectively (Fig. [3]), and Q p denotes 



In this section, we shall discuss our results for 
the chemical enrichment during the assembly of 
atomic-cooling halos and compare them with re- 
cent observations. We would like to emphasize 
that not all stars belonging to the Milky Way halo 
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Fig. 4. Predicted distributions (shaded areas) of low-mass 

stars in the [Fc/Ca] — [Ca/H] plane for /3 — 0.1. Upper panel: 
The full distribution. Filled (red) circles are observations of 
giants in the Galactic halo by Cayrcl et al. (2004), filled (red) 
squares are observations of dwarfs by Cohen ct al. (2004), filled 
(red) diamonds are three field stars depleted in n-capturc ele- 
ments (Ivans ct al. 2003), while open (red) triangles arc obser- 
vations of giants in dwarf spheroidal galaxies (Shctronc ct al. 
2001), and the (red) star is Draco 119 (Fulbright ct al. 2004). 
The dashed, thick and thin (blue) lines indicate, respectively, 
the 1(7 (innermost 68.3%) and 3cr (innermost 99.7%) cosmic 
scatter in the simulation. Lower panel: The partial distribu- 
tion of simulated stars for which > 90% of the total atmo- 
spheric Ca abundance is synthesized in PISNc. This distribu- 
tion is multiplied by a factor of 100. The observations and 
the dispersion of the full distribution arc plotted for compar- 
ison. Note, the mildly enhanced density of simulated stars at 
[Fc/Ca] ~ —1.6 originates from a nearly constant Fc/Ca yield 
ratio in the mass range 160 < m/ A4q < 170. 

were formed in atomic cooling halos, in particu- 
lar the stars at the metal-rich end (cf. Robertson 
et al. 2005). However, a significant fraction of 
the stars below [Ca/H] = —2 and a majority of 
the second generation stars presumably originated 
from this type of halo0 If so, they are likely to be 
found over the entire Galactic halo (Scannapieco 
et al. 2006). With this in mind, we will make no 
distinction in origin between our simulated stars 
and the field stars observed in the Galactic halo. 
Instead of iron, we choose calcium as the refer- 



![A/B] = log 10 {N A /N B ) t - \og 10 {N A /N B ) Q , where N x is 
the number density of element X. 



ence element for the metallicity of the stars. We 
do so for two reasons. First, the strong Ca II K res- 
onance line is commonly used as metallicity indi- 
cator (Beers et al. 1990) in surveys of metal-poor 
stars, such as the HK survey (Beers et al. 1992) 
and the Hamburg/ESO survey (Wisotzki et al. 
2000) . The much weaker Fe lines will remain unde- 
tected in low-quality survey spectra, although they 
can be detected by taking high-resolution spectra. 
This is of importance since the vast majority of 
the most metal-poor stars known thus were origi- 
nally recognized as being metal-poor on the basis 
of their Ca abundance, not their Fe abundance. 
The second reason we choose calcium is that the 
iron yield is, to a larger extent than calcium, af- 
fected by the unknown amount of fall-back in core 
collapse SNe and it varies by more than 3 orders 
of magnitude over the PISN mass range. The Ca 
yield, on the other hand, varies by merely 1 or- 
der of magnitude (Heger & Woosley 2002; Umeda 
& Nomoto 2002). For the employed observational 
data, either existing Ca abundances were used or, 
in case of the observed Galactic halo MDF, an off- 
set of [Ca/Fe] = +0.4 dex was applied to convert 
from Fe abundances to Ca abundances. This offset 
is based on the average Ca excess (Pagel & Taut- 
vaisiene 1995) observed in Galactic halo stars (see 
e.g., Cayrel et al. 2004), although there are out- 
lier stars in the lowest metallicity range which de- 
viate from the mean [Ca/Fe] = +0.4. Apart from 
the offset, however, the shape of the Galactic halo 
MDF probably doesn't change much whether Fe 
or Ca is used as the metallicity indicator. 

4.1. Comparison with Observed Metal- 
Poor Stars 

In Figure |4j we show the probability density 
function of low-mass stars in the [Fe/Ca] — [Ca/H] 
plane from a simulation in which 10% by num- 
ber (i.e. (3 = 0.1), or about 50% by mass, of all 
Pop III stars were allowed to explode as PISNe. 
The agreement between the model and the ob- 
servations of Galactic field stars by Cayrel et al. 
(2004) and Cohen et al. (2004) is satisfactory. It 
should be noted that the full distribution displayed 
in the upper panel of Figure 0] is effectively unal- 
tered for any f} < 0.1. 

The surprisingly small star-to-star scatter and 
the lack of evolution in the scatter that is ob- 
served in various abundance ratios for well de- 
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fined samples of extremely metal-poor Galactic 
halo stars (e.g., Cayrel et al. 2004; Arnone et 
al. 2005; Barklem et al. 2005) is not fully un- 
derstood. Such a small scatter could be an indi- 
cation of extremely short mixing timescales and a 
well mixed ISM already at very early times. How- 
ever, as evident from Figure |4j the small scat- 
ter is as well reproduced by our model over the 
entire mctallicity regime. Despite the relatively 
slow mixing assumed in our model, the small dis- 
persion is explained by a selection effect favoring 
contributions from SNe in a certain mass range 
(13 < mj M.q < 20, in this particular case) for 
the most metal-poor stars and by the averaging 
of a large number of contributing SNe at higher 
metallicities (Karlsson & Gustafsson 2005). Ul- 
timately depending on the stellar yields, a small 
star-to-star scatter can thus be realized without 
invoking unrealistically short mixing timescales in 
the early ISM. 

4.2. Characteristics of PISN-Dominated 
Stars 

In the lower panel of Figure 01 only those sim- 
ulated low-mass stars are displayed for which > 
90% of the total atmospheric calcium abundance 
was synthesized in PISNe. This regards the abun- 
dance originally present in the gas out of which the 
stars were formed and not the abundance at late- 
stage evolution, although the difference should be 
negligible for Ca. Henceforth, such stars will be 
referred to as PISN-dominated stars. From the 
results depicted in Figure HI three conclusions re- 
garding the chemical legacy of PISNe can immedi- 
ately be drawn. First and foremost, the low-mass 
stars with a dominant contribution from PISNe, 
as reflected by the Ca abundance, are located at 
significantly higher [Ca/H] than the most Ca de- 
ficient stars in the simulation, which are found at 
[Ca/H] ~ -4. The PISN-dominated stars are in- 
stead found in the range —3 < [Ca/H] < —2, 
with a maximum around [Ca/H] ~ —2.3. This 
prediction is at variance with the classical picture 
in which the chemical signatures of the first, very 
massive stars are anticipated to be found in the 
most metal-poor regime. This is exactly the effect 
illustrated in Figure [T] and described in ££2] Look- 
ing for the stars with the lowest abundances of Ca 
(or Fe) may thus not be the best way to find the 
relics of the first PISNe. 




-4.0 -3.5 -3.0 -2.5 -2.0 



[Ca/H] 

Fig. 5. — Predicted metallicity distribution functions 
as measured by [Ca/H] for f3 = 0.1 (dark blue lines) 
and /3 = 0.01 (light blue lines). Upper panel: The 
solid lines denote the total MDFs while the dashed 
lines denote the MDFs of PISN-dominated stars, for 
which > 90% of the Ca originate from PISNe. The 
binned data (red line) are the observed (not complete- 
ness corrected) distribution of Galactic halo stars from 
the Hamburg/ESO survey by Beers et al. (2005), 
where the vertical lines indicate the (la) stochastic 
uncertainty. Note that these data are shifted +0.4 
dex. Lower panel: The fraction of PISN-dominated 
stars. The majority of these stars fall above the upper 
limit of [Ca/H] = — 2.6 (dashed brown line), the range 
where the Hamburg/ESO survey is incomplete. 

Second, the fraction of PISN-dominated stars 
in the halo is predicted to be small. As illustrated 
in the lower panel of Figure 21 the distribution of 
PISN-dominated stars has to be multiplied with a 
factor of 100 in order to be detected on the same 
scale as the total distribution, shown in the up- 
per panel. Even for a relatively large (3 of 0.1 
(i.e., ~ 50% by mass), the total fraction of PISN- 
dominated stars below [Ca/H] = —2 is merely 
~ 2 x 10~ 3 , or about one star in 500. For lower (3, 
this fraction decreases further (see i j4.3D . 

Third, as noted in the chemical signature of 
PISNe distinctly differs from that of the less mas- 
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sive core-collapse SNe. This provides an opportu- 
nity to identify outliers, at relatively high metal- 
licities, as possible PISN-enriched, second gener- 
ation stars. As shown in Figured] (lower panel), 
PISN-dominated stars are predicted to be found 
both above and below the IMF-averaged value of 
[Fe/Ca] ~ —0.4. Noticeably, none of the stars in 
the samples by Cayrel et al. (2004) and Cohen et 
al. (2004) show a large enough deviation from 
the mean [Fe/Ca], and other abundance ratios, 
to be classified as potential PISN-enriched stars. 
However, the four stars belonging to the Draco, 
Ursa Minor, and Sextans dwarf spheroidal (dSph) 
galaxies have a relatively high [Fe/Ca] ~ —0.1. 
These gas-poor satellite galaxies primarily consist 
of old and metal-poor stars and were presumably 
formed at an early stage, as a result of the collapse 
of dark-matter halos of < 10 9 Mq. 

Could these dSph stars be second generation, 
PISN-dominated stars? This seems unlikely. The 
low a- element /iron ratios are generally attributed 
to Type la SNe in connection to lower SFRs in 
these systems (e.g., Matteucci 2003), and, while 
no contribution of s-process elements by asymp- 
totic giant branch stars is detected in dSph stars 
below [Fe/H] ~ -1.8 the high Ba/Y and low Y/Eu 
ratios (as well as low a /Ye) could be explained 
by the absence of elements synthesized in the a- 
process (Venn et al. 2004 and references therein), 
suggested to occur in hypernovae (Nakamura et al. 
2001). A few other, curious outliers such as Draco 
119 (Fulbright et al. 2004) and BD + 80° 245, 
G4-36, and CS 22966-043 studied by Ivans et al. 
(2003), will be further discussed in [J5] 

In the simulations, the distribution of stars with 
a dominant contribution from PISNe peaks around 
[Ca/H] ~ —2.5, or even higher (see Fig. [4] lower 
panel and Fig. O upper panel). In fact, the ma- 
jority of PISN-dominated stars are predicted to 
be found above [Ca/H] = —2.6 (corresponding to 
[Fe/H] = —3). This is the upper limit for which 
the Hamburg/ESO survey is considered to be com- 
plete (see Fig. [5l upper panel). A significant 
fraction of the second generation stars may there- 
fore have a chance to remain undetected in the 
Hamburg/ESO survey, not because they are too 
Ca-poor but because they are too Ca-rich. Note 
the deficit above [Ca/H] = —2.6 in the observed 
MDF which, to a large extent, has to be attributed 
to this intended observational bias against metal- 
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Fig. 6. — The predicted integrated (total) fraction of 
PISN-dominated stars below [Ca/H] = —2 as a func- 
tion of /3, corresponding to > 90% (big blue dots), 
> 99% (medium sized, blue dots), and > 99.9% (small, 
dark blue dots) PISN-enrichment. The dashed (red) 
lines indicate the observational upper limit of f3, as- 
suming that none of the ~ 600 Galactic halo stars with 
[Ca/H] < —2 for which high-resolution spectroscopy is 
available show any signature of PISNe (N. Christlieb, 
priv. comm.). The dotted (blue) line and shaded 
(blue) area denote the predicted range of antici- 
pated from the data of Padoan et al. (2007) and Greif 
& Bromm (2006), respectively. If nothing else stated, 
it is assumed that the PISN signature is detectable in 
Galactic halo stars with > 90% PISN-enrichment. 

rich stars (Fig. upper panel). With this taken 
into account, the observed and predicted MDFs, in 
particular the one with (3 — 0.01, agree to within 
the stated uncertainties. The lower panel of Fig- 
ure [5] shows the fraction of PISN-dominated stars 
for = 0.1 (dark blue line) and f3 = 0.01 (light 
blue line). Evidently, the fraction also peaks sig- 
nificantly above the cut-off of the observed MDF 
at [Ca/H] ~ —4 and goes to zero approaching 
this cut-off. We have additionally run a simula- 
tion with iron as the reference element. The cor- 
responding distribution of PISN-dominated stars 
(here defined such that > 90% of the Fe is syn- 
thesized in PISNe) peaks at [Fe/H] ~ —2.5, which 
concurs with the above result that these stars are 
to be found at relatively high metallicities. In this 
simulation, there is a weak tail of PISN-dominated 
stars towards low metallicities, originating from 
the least massive PISNe. However, for (3 = 0.01, 
merely ~ 1 star in 1000 show a clear PISN signa- 
ture at [Fe/H] = —4, nearly a factor of 100 less 
than around the peak at [Fe/H] = —2.5. 
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4.3. Observational Constraints on Pop III 
Supernovae 

As it appears from Figure the fraction of 
PISN-dominated stars depends strongly on the pa- 
rameter j3, the fraction of Pop III stars becoming 
PISNc, which is to be expected. This dependence 
may be used to put constraints on the number of 
PISN-dominated stars expected to be found in the 
metal-poor Galactic halo and, conversely, to put 
observational upper limits on (3 itself. Figure [6] 
shows the integrated fraction of PISN-dominated 
stars as a function of [3. This fraction, termed a^, 
is defined as the number of PISN-dominated stars 
below [Ca/H] = — 2 divided by the total number of 
stars below [Ca/H] = —2. The fractions displayed 
in Figure [6] are based on the assumption that a 
distinct PISN signature is present over the entire 
PISN mass range. In the reminder of the paper, 
we shall further explicitly assume that this signa- 
ture is detectable in stars if, and only if, > 90% 
of the atmospheric Ca abundance originates from 
the PISNe. In Figure[6l we have also indicated the 
predicted fractions corresponding to 99% and 
99.9% PISN-enrichment (see also Table [I]). These 
results will be briefly discussed later in this sec- 
tion. 

The number of observed Galactic halo stars 
below [Ca/H] = —2 for which abundance analy- 
ses based on high-resolution spectroscopy exists is 
roughly 600 (N. Christlieb, priv. comm.). None 
of these ~ 600 stars appear to show a dominant 
PISN signature. If so, we estimate that (3 < 0.07 
(see Fig. |6|), which corresponds to a mass frac- 
tion of very massive Pop III stars of < 40%. Here, 
we disregard Draco 119 (Fulbright et al. 2004) 
and the group of field stars studied by Ivans et 
al. (2003) as PISN-dominated, second generation 
stars (see ij5|). An upper limit of j3 = 0.07 is con- 
sistent with the predicted range of derived by 
Greif & Bromm (2006) and the j3 derived from 
the IMF by Padoan et al. (2007), which all fall 
in the range 0.005 < (3 < 0.02. Conversely, as 
indicated in Figure El this range corresponds to a 
range in of 1.5 x 10 -4 < < 5 x 10~ 4 , which, 
at face value, means that ~ 3 — 10 times as many 
metal-poor (i.e., [Ca/H] < —2) Galactic halo stars 
need to be analyzed in order to find a single second 
generation star with > 90% of its Ca originating 
from primordial PISNe. In these results will 
be discussed further. 



Table 1: Fraction of PISN-enriched stars as a func- 
tion of level of PISN-cnrichmcnt 



Level of 


Fraction of 


PISN-enrichment a 


PISN-enriched stars b 


> 1% 


6.97 X 10~ 2 


> 10% 


2.34 x 10~ 2 


> 50% 


3.60 x 10~ 3 


> 90% 


3.04 x 10~ 4 


> 99% 


7.35 x 10" 5 


> 99.9% 


6.75 x 10~ 5 







a As measured by the amount of Ca originating from PISNe. 
b All fractions are calculated for /3 = 0.01. 

Also shown in Figure [6] are the predicted frac- 
tions of stars that are pre-enriched to at least 99% 
(medium sized, blue dots) and 99.9% (small, dark 
blue dots) by PISNe. These will be referred to as 
purely PISN-enriched stars. Although the number 
of such stars is smaller than the number of stars in 
which up to 10% of the metals (Ca) are allowed to 
originate from core collapse SNe, it is not radically 
smaller. In the range 0.005 < (3 < 0.02, e.g., the 
predicted fraction of purely PISN-enriched stars is 
only smaller by factor of ~ 4 (see also Table [1} . 
Moreover, assuming that none of the ~ 600 stars 
below [Ca/H] = —2, for which high-resolution 
spectroscopic data are available, are purely PISN- 
enriched stars, the range of possible values for (3 
is constrained to be < 0.2. As indicated above, 
a slightly stronger constraint may be placed on (3 
if the reasonable assumption is being made that 
none of these stars has a PISN-enrichment which 
exceeds 90%. However, due to the increasing risk 
of making false identifications, [3 may not be con- 
strained much further by allowing for lower levels 
of PISN-enrichment than 90%. Table [Q shows the 
predicted fraction of stars formed from gas pre- 
enriched by PISNe to various levels, for (3 — 0.01. 
Clearly, the fraction of PISN-enriched stars in- 
creases with decreasing level of PISN-enrichment. 
However, even though the expected number of 
stars with, e.g., a PISN-enrichment of > 50% is 
~ 10 times more than the corresponding num- 
ber with > 90% PISN-enrichment, we would in 
practice not gain much as it would most likely 
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be very difficult to correctly identify many of 
the less PISN-dominated stars. Note also that 
for (3 = 0.01, less than 10% of all stars below 
[Ca/H] = —2 are formed from gas enriched to any 
level by PISNe. 

4.4. Parameter Dependence and Sensitiv- 

ity 

We have performed a number of simulations 
where we varied the values of several model pa- 
rameters, in order to illustrate the sensitivity of 
our results to the uncertainties in these parame- 
ters. The SFR is a crucial ingredient in the model 
and we have varied both the Pop III and Pop II 
SN rates by changing un,o and un,0j respectively. 
The turbulent diffusion coefficient D t in the early 
ISM is difficult to estimate. According to Bateman 
& Larson (1993), the diffusive mixing of matter 
in the cold, neutral medium in the Galactic disk 
(e.g., through molecular cloud collisions) occurs on 
a similar time scale as the one adopted here, while 
the mixing should be substantially faster in hot, 
ionized media. Alternatively, Pan & Scalo (2007) 
developed a "slow" mixing model with inefficient 
(i.e., low D t ) diffusion to explain the apparent ex- 
istence of primordial gas at relatively low (z ~ 3) 
redshift. We have varied D t in our model to map 
out the effect. 

The results of the simulations are summarized 
in Table [2 which shows the dependence of a^y on 
um,o, 1*11,0, and D tl respectively. In all simulations, 
(3 is fixed to 0.01. The impact on a^y when altering 
the Pop III (col. 1) and Pop II (col. 2) SN rates 
is non-trivial, with, respectively, a local maximum 
and a local minimum in the value of appear- 
ing within the examined ranges of the values for 
these rates. The dependence on D t (col. 3) is less 
complex. We shall not discuss these dependences 
in detail but notice that a-yy is roughly propor- 
tional to J Mm undt / / u\dt and that the Pop III 
and Pop II rates are coupled through Q p (see eqs. 
©-©)■ In particular, the counter-intuitive de- 
crease in Oyy with an increasing Pop III SN rate 
above u M ,o = 1-33 x 10~ 2 kpc~ 3 Myr" 1 (= 1/3 
of the adopted um,o) may be understood from the 
fact that an initially high Pop III rate more effi- 
ciently suppresses subsequent formation of Pop III 
stars, while it simultaneously promotes the forma- 
tion of Pop II stars, as a consequence of the change 
in the fraction of primordial gas. Note that a sim- 



Table 2: The dependence of a-yy on model parameters. 
The three parameters which were varied are listed in the 
top row; from left to right, these are the Pop III SN rate, 
the Pop II SN rate, and the turbulent diffusion coefficient. 
The leftmost column shows the factor by which these pa- 
rameters were changed from their fiducial values, and the 
resulting values for a T/ are shown in the three rightmost 
columns below the specific parameter which was varied. 



Parameter 3, 


«M,0 




«n,o 




Dt 




X10 


1.07 x 10" 


-4 


7.35 x 10" 


-4 


6.89 x 10- 


5 


X 3 


1.57 x 10- 


-4 


1.99 x 10" 


-4 


7.98 x 10" 


5 


X l b ' c 


3.04 x 10" 


-4 


3.04 x 10" 


-4 


3.04 x 10" 


4 


x 1/3 


4.78 x 10" 


-4 


4.59 x 10" 


-4 


1.02 x 10" 


3 


x 1/10 


4.01 x 10" 


-4 


3.41 x 10" 


-4 


1.92 x 10" 


3 







a The fiducial values of the model parameters are discussed in 

tf51 All fractions are calculated for fS = 0.01. 
b a- rl = 4.78 X 10~ 4 , adopting the set of PISN yields calculated 

by Umeda & Nomoto (2002). 
c a-fy = 6.03 X 10 -4 , assuming that SFR <x p 2 



ilar effect occurs when j3 is increased, indicated in 
Figure [3] 

The overall change in over the entire two 
decade range in un,o, as well as in un,0i is merely 
a factor of ~ 4, which is quite moderate. The 
sensitivity on D t is higher, a factor of ~ 30 over 
the examined range. We also ran one simulation 
where the SFR was set proportional to the square 
of the density instead of the normal linear depen- 
dence, and one where we changed the set of PISN 
yields to those by Umeda & Nomoto (2002). In 
each of these two runs, a-yy changed by less than a 
factor of 2 (see footnote b and c in Table ^ . 

The data presented in Table [2] are a measure of 
the uncertainty in our calculations. Evidently, the 
result is insensitive to the SFR and PISN yields, 
while it is relatively sensitive to the amount of 
mixing in the ISM. If we, e.g., adopt a 10 times 
smaller turbulent diffusion coefficient, the fraction 
of PISN-dominated stars is increased by a factor 
of 6. Still, even for such a low D t , the range of 
f} predicted by Greif & Bromm (2006) would be 
consistent with the observational upper limit. We 
would like to emphasize that in all the runs that 
we have carried out, the vast majority of the PISN- 
dominated stars have significantly higher Ca abun- 
dances than the most Ca-poor stars in the simu- 
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lations, irrespective of the predicted value of a™. 
Moreover, for all the parameters that we have var- 
ied, (3 remains the parameter to which the results 
are most sensitive (cf. Fig. The above findings 
indicate two points: first, due to the sensitivity on 
(3, we have the means to measure the frequency of 
very massive stars in the early universe from ob- 
servations of large, homogeneous samples of metal- 
poor stars in the Milky Way halo. Second, due to 
the insensitivity on model parameters apart from 
/3, our results and the conclusions based on these 
results should be relatively robust. 

5. Discussion and Summary 

We have shown that Milky Way halo stars born 
in atomic cooling halos of M ~ 10 8 Mq and dis- 
playing a chemical abundance pattern character- 
ized by that of Pop III PISNe are to be found at 
significantly higher metallicities than where the 
most metal-poor stars have been found. Conse- 
quently, although stars in the metal-poor tail of 
the MDF, like the extremely iron-deficient stars 
HE 0107-5240 (Christlieb et al. 2002) and HE 
1327-2326 (Frebel et al. 2005), are very important 
for the understanding of the star formation pro- 
cess (e.g., Bromm & Loeb 2003; Umeda & Nomoto 
2003; Tumlinson 2007) and feed-back mechanisms 
(e.g., Karlsson 2006) in the early universe, they 
seem to tell us little about primordial PISNe. The 
majority of Pop II stars with a dominant contri- 
bution from PISNe are predicted to have Ca abun- 
dances in excess of [Ca/H] = —2.6, the upper limit 
for which the Hamburg/ESO survey of metal-poor 
stars is considered to be complete. This implies 
that a significant fraction of the PISN-dominated 
stars may have escaped detection, not because 
they are too Ca-poor but because they are too 
Ca-rich. The number fraction of very massive 
Pop III stars exploding as PISNe is estimated to 
be (3 < 0.07 (< 40% by mass), assuming that cur- 
rently no Galactic halo star below [Ca/H] = — 2 
with available detailed abundance analysis shows 
a chemical signature characteristic of PISNe. This 
result is consistent with the theoretical prediction 
of 0.005 < (3 < 0.02, which corresponds to a hand- 
ful of PISN-dominated, second generation stars in 
a sample of 10 000. The number of purely PISN- 
enriched stars is predicted to be a factor of ~ 4 
smaller, assuming this range of (3. 



In an interesting paper, Salvadori et al. (2007) 
also discuss the existence of second generation 
stars and their location in the MDF. Their 
stochastic model of chemical enrichment differs 
from ours in a number of aspects. Firstly, they 
model the chemical evolution of the Galaxy as 
a whole based on the merger-tree approach (cf. 
Tumlinson 2006), while we focus on the evolution 
of atomic cooling halos using the mixing volume 
picture. They also assume that all Pop III stars 
are very massive, while we allow for less mas- 
sive Pop III star formation as well, controlled by 
the parameter (3. Hence, in their models with 
^crit > (for Z cr it = 0, no PISNe are formed), all 
second generation stars are formed from gas only 
enriched by PISNe. Finally, we note that while 
our chemical enrichment model accommodates the 
important contribution of metals ejected by the 
supernovae from Pop III stars formed in miniha- 
los, with masses less than 10 8 Mq (see Bromm et 
al. 2003; Greif et al. 2007), this contribution is 
neglected in Salvadori et al. (2007) based on the 
assumption that strong photodissociating radia- 
tion from the first stars quenches their formation. 
It appears likely, however, that this effect may 
not have substantially lowered the Pop III star 
formation rate (e.g. Ricotti et al. 2001; John- 
son et al. 2007; O'Shea & Norman 2007; Wise & 
Abel 2007; see also Yoshida et al. 2003). Despite 
these differences, both models generate broadly 
similar MDFs. However, at variance with the re- 
sults presented here, Salvadori et al. predict that 
the highest fraction of PISN-dominated stars is 
to be found in the lowest metallicity regime and 
only a very small fraction at [Fe/H] > —1, pro- 
vided that the host halo is able to form both first 
and second generation stars. This finding may 
be a consequence of these authors' assumption of 
the instantaneous mixing of metals in the galactic 
medium, but should be analyzed in more detail to 
be fully understood. 

The most characteristic and reliable chemical 
signature of PISNe is probably the absence of 
neutron-capture elements. Although this may not 
be a unique fingerprint (massive core collapse SN 
may be unable to synthesis n-capture elements as 
well) , we have shown that stars formed from mate- 
rial pre-enriched by PISNe emerge in a metallicity 
regime where both the r-process and the s-process 
signatures are thoroughly established (see Burris 
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et al. 2000). The presence of n-capture elements 
in normal Pop II stars will thus serve to rule out 
PISN pre-enrichment. Metal-poor stars with en- 
hanced abundances of r-process elements are cer- 
tainly interesting in their own right as they, e.g., 
can be used as cosmochronometers and allow us 
to study the origin and nature of the r-process. 
However, as important as it is to search for r- 
process enhanced stars (see e.g., Christlieb et al. 
2004), alternatively, the search for stars unusu- 
ally depleted in r-process elements may help us 
locate candidates of the "missing" second stellar 
generation. With high-resolution follow-up spec- 
troscopy, further signatures indicative of PISNe 
such as a pronounced odd-even effect and, e.g., 
an atypical Fe/Ca ratio (see Fig. [4]), should then 
be measured to ultimately verify or refute candi- 
dates as PISN-dominated, second generation stars. 
Although the actual number of PISN-dominated 
stars is predicted to be small, perhaps only one 
or two in 10 000 below [Ca/H] = —2, the ongoing 
and future large-scale surveys, such as the SDSSq 
SEGUEi, and LAMOST0, aim to find at least an 
order of magnitude more stars in this metallicity 
regime than is presently known. A few thousand 
of these stars should be bright enough to be suc- 
cessfully observed with high-resolution by 8 — 10 m 
class telescopes (Beers & Christlieb 2005). With 
the advent of 30 — 100 m class telescopes, high- 
resolution spectroscopic follow-up will become fea- 
sible for the majority of stars discovered below 
[Ca/H] = -2. 

There are a few examples of metal-poor stars in 
the Galactic halo which do show anomalous abun- 
dance patterns. Ivans et al. (2003) studied a group 
of three dwarf/subgiant stars with [Ca/H] ~ —2.2, 
BD 80° 245, G 4-36, and CS 22966-043, which 
appears to be heavily depleted in n-capture ele- 
ments (see also, e.g., HD 122563, recently studied 
by Mashonkina et al. 2008). These stars also dis- 
play low abundances of Na, Al, and a-elements, 
as indicated by the relatively high [Fe/Ca] ratios 
in Figure H Interestingly, BD 80° 245 and G 4- 
36 both show a particularly low [Al/Ca] ratio and 
BD 80° 245 and CS 22966-043 have [Mg/Ca] < 0. 



2 The Sloan Digital Sky Survey 

3 The Sloan Extension for Galactic Understanding and Evo- 
lution. Operating since 2005. 

4 The Large sky Area Multi-Object fiber Spectroscopic Tele- 
scope. First spectrum delivered May 28, 2007. 



However, the [Na/Ca] ratio in all three stars is 
found to be typical for Galactic halo stars and al- 
though the heavy iron-peak elements are low in 
BD 80° 245, they are strongly enhanced in G 4- 
36, and CS 22966-043, an observation which is not 
understood. Ivans et al. (2003), explain the pecu- 
liar abundance patterns as a result of an early con- 
tribution by thermonuclear SNe (Type la) in the 
subsystems in which these stars were born. They 
do not consider PISNe as an alternative, presum- 
ably because very massive stars traditionally are 
believed to contribute to the chemical enrichment 
only at the very lowest metallicities. As intrigu- 
ing as it is to envision these objects, in particular 
BD 80° 245, as PISN-enriched, second generation 
stars, the match is ambiguous. Further investiga- 
tion should, however, be conducted on these pecu- 
liar stars. Another star shown in Figure|4]is Draco 
119. This [Fe/H] = —2.9 giant is a member of the 
Draco dSph and is the only star currently known 
to be completely devoid of elements above Ni, to 
within detection limits (Fulbright et al. 2004). 
The upper limit on the Ba abundance, e.g., is esti- 
mated to [Ba/Fe] < —2.6, more than a full dex be- 
low the next most Ba-poor stars at [Fe/H] ~ —3. 
Based on the light element abundances, Fulbright 
et al. (2004) conclude that Draco 119 was born out 
of gas primarily enriched by massive core collapse 
SNe, and not by PISNe, which is consistent with 
our finding that PISN-dominated stars generally 
should appear above this metallicity regime. 

We have assumed that the metal-poor Galactic 
halo mainly consists of stars formed during the as- 
sembly of atomic cooling halos. We note, however, 
that a significant fraction of these halos may have 
ended up in the bulge of the Milky Way. With 
the new high-resolution infrared spectrographs on 
the 8 — 10 m class telescopes, detailed spectral 
analysis of dust-obscured Galactic bulge stars has 
recently become feasible (e.g., Ryde et al. 2007). 
The search for relics of the first stars should there- 
fore be pursued also in this, previously inaccessible 
part of our Galaxy. 
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